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Howard University , Washington, D.C . 20059 
and 

Thomas W. Kerslake§ 
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A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat 
receivers employing encapsulated phase change materials (PCMs) that undergo cyclic melting and freezing. The 
present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit 
(LEO). The solar heat receiver is integrated into a closed -cycle Brayton engine that produces electric power during 
sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both 
on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power 
crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to 
the working fluid are derived, all of which are related to the SOC through measurable parameters . Lower and upper 
bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, 
geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. 
Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis 
Research Center’s Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of 
measurable parameters. 


Nomenclature 

A = area or growth constant 

c = specific heat of solid or liquid PCM 

c p = specific heat of working fluid 

D cav = active cavity diameter 

D ap = aperture diameter 

F = geometric view factor 

h = enthalpy per unit mass 

h sf = PCM latent heat of fusion 
H = Heaviside function 

m = working fluid mass flow rate 

M = total number of axial nodes along tube 
or total PCM mass 

N = total number of tubes in receiver 

p = working fluid pressure 

Q = heat transfer rate 

R = gas constant 

s,S = specific, total entropy 

S gen = entropy generation rate 
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Ste = Stefan number 

t = time 

T = temperature 

T m = PCM melting temperature 

T 0 = environmental dead state temperature 

T p ,T t = sunset, sunrise temperature 

T* = virtual source temperature 

u,U = specific, total internal energy 

V = total volume 

W = rate of work transfer 

z = axial location 

Pj = first conjugate SOC function 

p 2 = second conjugate SOC function 

Xj = jth tube mass fraction 

£ = thermal capacitance ratio 

O = primary SOC function 

y = ratio of specific heats 

p = density 

a = Stefan-Boltzmann constant 

T on , T 0 ff = sun period, eclipse period 

T cyc = total orbit period 

Subscripts 

avg = average 

in, out = tube inlet, tube outlet 

losses = losses through shell and aperture 

min, max = minimum, maximum 

rcvr = receiver 
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Fig. 1 Thermodynamic cycle for closed Brayton engine integrated with solar heat receiver. 


Introduction 

S OLAR heat receivers are very critical components in 
the production of electric power via solar dynamic 
power systems (SDPSs). During operation, the SDPS 
uses: 1) a concentrator to collect and focus the incident 
energy onto the aperture plane of a central receiver, 2) a 
central receiver to collect and distribute, with minimal 
losses, the reflected energy from the concentrator, 3) 
working fluid tubes aligned along the periphery of the 
receiver to absorb the distributed energy as heat, thus, 
raising the temperature of the working fluid (typically a 
low-Prandtl-number fluid) flowing through the tubes, 4) 
a turbine to expand the high temperature working fluid 
to produce mechanical work via a rotating shaft, 5) a 
compressor to circulate the working fluid through the 
working fluid tubes, and 6) an alternator to convert 
mechanical shaft motion into electric power. A 
recuperator is often added to increase the thermal 
efficiency of the thermodynamic cycle (typically a 
closed Brayton cycle as depicted in Fig. 1). 

Solar heat receivers employing encapsulated phase 
change materials (PCMs) have the advantage over 
sensible heat receivers of requiring less mass while 
producing higher energy storage densities. This, in 
turn, makes them ideal candidates for energy storage in 
the space environment where temperatures are 
sufficiently high and PCMs with high latent heats of 
fusion become indispensable. 

In this paper, a theoretical framework on the so-called 
thermal state -of -charge (SOC) of solar heat receivers 
employing latent heat thermal energy storage (LHTES) 
is developed. The instantaneous amount of phase 


change material (PCM) in the liquid phase was 
identified by Strumpf et al. 1 as an indicator of the SOC. 
This definition, however, is a better indicator of PCM 
effectiveness or some performance measure (e.g. 
efficiency) of the receiver as it relates to incorporating 
phase change storage. It may also be tempting to define 
the SOC as the instantaneous amount of energy stored in 
the receiver. According to NASA 2 , “ Techniques are 
needed to determine the so-called receiver state of 
charge, or the quantity of stored thermal energy within 
the receiver .” However, this idea can be quickly 

dismissed on second law grounds since energy quality 
can be considered a factor in determining the true SOC. 
It should be pointed out that the issue of energy quality 
does not adhere to conservation principles. In other 
words, the statement conservation of entropy has no 
meaning since all real devices that undergo energy 
exchange processes are involved in the one-way 
production of entropy. In what follows, it will be 
shown that the available power stored in the receiver is 
related to a newly-defined, time-dependent SOC 
function, which may be completely characterized by 
measurable parameters. Knowledge of the SOC allows 
for better control strategies relating to power 
management schemes during such operations as peak 
power demand and emergency shutdowns with 
subsequent restarts. It also helps to better identify the 
energy startup characteristics of the solar heat receiver 
in relation to the entire solar dynamic (SD) system, 
which ensures safe operation of the SD system through 
all modes and regimes of operation. 
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Fig, 2 Solar heat receiver available power and SOC model indicating control volume. 


Theoretical Framework 

Solar Heat Receiver Available Power 

By definition, the available power of any device is the 
maximum rate at which energy may be extracted by a 
work transfer interaction if the device is allowed to 
come into total (thermal, mechanical, chemical) 
equilibrium with its surroundings at some dead state. 3 
Shown in Fig. 2 is the model (including control volume) 
used to derive an expression for the available power 
stored in the receiver. A 1 st law energy balance on the 
entire receiver with a single fluid stream results in the 
following: 

-it (1) 

where W is the rate of work transfer across the 
boundary of the control volume (this is what could be 
theoretically extracted if the receiver was connected to a 

work-extracting mechanism), m is the working fluid 
mass flow rate, h is the enthalpy per unit mass of the 

working fluid, Q rcvr is the rate at which energy crosses 


the aperture plane, Q hsses is the rate at which energy 

leaves the receiver due to reradiation from the canister 
surfaces back out through the aperture and conduction 
losses through the receiver shell, and U is the total 
internal energy of the receiver. An associated entropy 
balance results in 

^ ^ ‘ ^ Qrcvr Q losses , o (^)\ 

a7 = rl-r s L + ~~ + ^ ( } 

where S is the total entropy of the receiver, s is the 
entropy per unit mass of the working fluid, T* is a 
virtual source or effective aperture temperature (defined 
in the next section), T 0 is the environmental dead state 

temperature, and S gen is the rate of entropy generation 
inside the receiver. Subsequently eliminating the power 
loss term between Eqs. 1 and 2 gives 


w = 

m (h~T 0 s) 

— 

m{h-T 0 s) 

+ Q 
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-T o S gen -j{U~T o S) (3) 


NASA/TM— 1 998-207920 


3 



in which the maximum is 



m(h-T a s) 

- 

m(h-T 0 s) 

+ G 

revr 

( J' > 

1 T 

T 


- 

in 

~ 

out 

V 1 J 


-^(U-T„S) W 


since S gen =0 for a receiver operating reversibly. 

Now, it is assumed that the specific enthalpy in Eq. 4 is 
a function of temperature and pressure, i.e. 


i V . N 

_ m, suc h that m = Y mi and V £ = 1 (12) 
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Eq. 1 1 can be expressed in non-dimensional form as 
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h = h(T,p) (5) 

and the specific entropy is a function of specific 
enthalpy and pressure, i.e. 

s = s(h,p) (6) 

which for changes in specific enthalpy and specific 
entropy results in 


dh = 


dh_ 

dT 


op 


dp 


(7) 


ds = 


05 

dh 


dh + 


05 

dp 


dp 


( 8 ) 


Through the use of Maxwell’s relations, ideal gas 
assumptions for the working fluid, and the definition of 
specific heat at constant pressure, it can be shown that 
Eqs. 7 and 8 when integrated from inlet conditions to 
outlet conditions yield 


Ku, ~K = C p( T ou, ~ T in) 


S ou,- S in = C p ln l 


/ T A 

1 out 

T 

V m J 
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r p 

r out 


V P,n J 


(9) 

(10) 


in which c p is the working fluid specific heat at constant 
pressure, y is the ratio of specific heats (Cp/c v ), P out is 
outlet pressure, and P in is inlet pressure. Furthermore, 
the internal energy U and entropy S are given, 
respectively, by 

v=±± !Sj(H, d K <> 4 > 

7=1 -= I V, 

and 

s-t± JIM,". <l5) 

J =l i=] Vg 

in which the integration takes place over the ith region 
and jth tube. Upon further defining the dimensionless 
parameters 
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where T = X on + T off , the following dimensionless 
receiver available power results: 


For a solar heat receiver with N tubes (see Fig. 2), the 
available power is written as 
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where upon defining the jth tube mass fraction as 
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where Ste = cT fh f I s the Stefan number, which is the 
ratio of PCM sensible heat to latent heat, and the ratio 
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Mc/^mt c j is t ^ erma ^ capacitance ratio expressing 

the relative amounts of sensible heat capacity of the 
PCM to sensible heat capacity of the working fluid. In 
addition, the dimensionless available power equation 
can be interpreted physically as follows: the first 

bracketed term is the available power loss from the 
receiver to the gas; the second term is the available 
power gain by the receiver due to the net heat 
interaction across the aperture plane; the last term is the 
available power loss or gain due to unsteady charging 
and discharging. Furthermore, an interesting 
comparison can be made between the fraction of 
incident power available at the aperture plane, 
expressed in Eq. 17 as 


T 1 

y/ = \ °- = \ 

T Y * t ** 


(18) 


and that which is reported in Moynihan 4 , i.e. 


4 T 1 (T } 
^’=1--— + - — 

3 T* 3^7* J 


=i-i- 1 


tf i 


3 T** 3 [T** 


(19) 


A more detailed graphic comparison is shown in Fig. 3. 
It should be pointed out that over the anticipated region 
of operation (0.1 < T</T* < 0.3) the agreement is good. 



Dead State to Virtual Source Temperature Ratio 


Fig. 3 Comparison of the fraction of incident power 
available at the aperture as a function of the ratio of 
dead state temperature and virtual source temperature. 


Virtual Source Temperature 

The effective temperature of the aperture due to the 
net heat interaction across the aperture plane is defined 
as the virtual source temperature. It expresses 
continuity of energy reradiated from the outer surfaces 


of the canisters to the aperture. Therefore, an energy 
balance on the aperture plane of the receiver shown in 
Fig. 2 results in 

s <o - r ' 4 w] = (o - t :] (20) 

;=1 

where upon solving for the virtual source temperature 
T* gives 


T (f) = 




;=i 
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( 21 ) 


which in dimensionless form is written as 
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In Eq. 22, the jth area ratio (Aj/A ap ) can be written as 



TtD A z 





V 


A z 


A 

) 


(23) 


Notice that Eq. 23 contains one of the cavity aspect 
ratios (D^y/Dap), which is a key parameter that affects 
the thermal performance of the solar heat receiver. 
Also, the geometric view factors in Eqs. 20-22 are given 
by analytic expressions found in Howell. 7 


Gas Available Power 

Recall the expression given by Eq. 11, which is the 
instantaneous available power stored in the receiver. It 
should be pointed out here that the first term in brackets 
in Eq. 1 1 represents the instantaneous available power 
of the gas before mixing in the outlet manifold. This 
power, which is the difference between the enthalpy 
transferred to gas and a term proportional to the entropy 
transferred to the gas, is rewritten here as 


W gas — C p 


(T„-T„\-TJn -S- 


+r 


(r - 1 
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V ^ „ 


( 24 ) 


which is expressed in dimensionless form as 
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which can be further written in the following non- 
dimensional form: 


Physically, the reduction in available power associated 
with the transport of entropy to the gas can be attributed 
to two sources: 1) heat transfer to the gas across finite 
temperature differences and 2) frictional effects leading 
to reductions in pressure along the lengths of each of 
the tubes in the receiver. The concept of entropy is 
associated with the amount of unavailable energy within 
a system. Therefore, the available power of the gas, 
expressed by Eq. 25, is that which is delivered to the 
oudet manifold before any mixing takes place. Any 
further reduction in available power takes place in the 
outlet manifold due to irreversible mixing of each of the 
individual fluid streams. Finally, it is observed that all 
of the parameters in Eq. 25 can either be readily 
calculated or directly measured. 

Mixing-Based Lost Available Power 

It is known that the lost available power associated 
with the gas is proportional to the entropy generation 
rate, where the proportionality constant is the dead state 
temperature. 3,6 Expressed mathematically, 

W lost =T 0 S g en (26) 

For a solar heat receiver with N tubes and N associated 
fluid streams, the entropy generation rate due to 
irreversible fluid stream mixing in the outlet manifold is 
given by 


n 

S gen = m s N+l (27) 

7=1 

which, due partly to Eq. 10, can be subsequently written 
in terms of temperatures as 


fr 


s gen =mC p '£ l Xj ln 

/= I 


(28) 
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Wu>« 

mc p T a 


7=1 


f T ^ 
*N ± I 


(30) 


As expected, when each tube in the receiver is imparted 
with the same incident flux, the temperature of each gas 
stream exiting all the tubes is the same, resulting in no 
loss in available power. This is revealed in Eq. 30. 

Relation Between SOC and Available Power 

A dimensionless conjugate SOC function [Jj is 
defined here as the ratio of instantaneous available 
power stored in the receiver with no available power 
lost to the working fluid to minimum gas available 
power required to operate the turbine, or 
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A second dimensionless conjugate SOC function p 2 is 
defined as the ratio of instantaneous gas available 
power to minimum gas available power required to 
operate the turbine, i.e. 


txy 
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where the outlet manifold mixing losses have been 
neglected in both conjugate functions for convenience. 
Upon further defining the denominator of Eqs. 31 and 
32 as 


where T } is the fluid outlet temperature of the jth tube 
just before entering the outlet manifold and N+l 
corresponds to mixed mean properties in the outlet 
manifold just before entering the turbine. The 
corresponding loss in available power of the gas is 
expressed as 


=i>, 

7=1 


(t* 
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V 1 


An A 
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(33) 


the conjugate SOC functions and P, mn can be related to 
the dimensionless receiver available power (Eq. 1 7) as 
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Pmi„(Pl— P2 ) = ^T (34) 

mCpT 0 

Now, define the primary SOC function as the 
dimensionless combination 


Q> = 


A-1 


(35) 


such that <I> is always in the range 0 < O < 1 . Notice 
that pmin is p 2 evaluated at T out = and P out = and 
Pmax is the maximum value that the first conjugate SOC 
function Pi can take on, which can be shown to be 


0- 



f v 
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c„ Ste 



txMTL-Kl-l n 


T L 



(36) 


mode (BOM) wherein measured quantities are 
repeatable (within allowable limits) from one sunrise to 
the next and from one sunset to the next. In order to 
describe the aforementioned minima in BOM, the gas 
inlet temperature profile must be specified. Owing to 
the cyclic nature of the solar source, thermodynamic 
parameters in the system responding to the cyclic solar 
source will also experience cyclic changes throughout 
the orbit cycles. It is interesting to note that the only 
coupling that the receiver has with the rest of the 
thermodynamic cycle is through the gas mass flow rate 
and gas inlet temperature, which is intricately coupled 
to the gas mass flow rate and components of the heat 
rejection loop (recuperator, gas coolers, radiators, heat- 
rejection coolant, etc.). Knowledge of these two 
measurable parameters along with the pressure drop 
through the tubes are all that is required to calculate the 
gas available power (Eq. 24). In order to model cyclic 
changes that occur inside the receiver in BOM, several 
inlet temperature test profiles are proposed: 

Sawtooth: 


which is just the steady-state equivalent of Eq. 31. 
Therefore, the maximum SO C is achieved when the 
thermodynamic state of the receiver is driven to steady- 
state conditions even though the system is designed to 
operate under cyclic conditions. 


Results and Discussion 

The intrinsic coupling between the size (and design) 
of the solar heat receiver and the turbine for which it is 
intended to supply high temperature, high pressure gas 
makes it prudent to understand the minimum necessary 
thermodynamic requirements for operating the turbine. 
Mason 8 describes a process called motoring in which a 
DC electric power source is initially used to drive the 
turbo-alternator compressor (TAC) while the turbine is 
pre-heated and, ultimately, becomes self-sustaining. 
Mason 8 identified the cycle temperature ratio (turbine 
inlet/compressor inlet temperature) as the leading 
indicator for the time when the TAC becomes self- 
sustaining. It was determined that minimum motoring 
time is achieved as the cycle temperature ratio 
approaches a value of three (3) asymptotically. This 
information can subsequently be used to determine the 
minimum thermodynamic state-point (temperature and 
pressure) and associated minimum gas available power 
needed to operate the turbine. Of course, this minimum 
gas available power is that which is delivered from the 
receiver to the gas. 

One important mode of operation of the solar 
dynamic (SD) system is the so-called balanced orbit 


= , 0</<T„ 

^ on 

T 0ff 1 off 


which can be written in dimensionless form as 



where T 



T +T 


off 


, Tp is the sunset 


temperature, and Tj is the sunrise temperature; T p and T] 
are repeatable from one cycle to the next in BOM. 
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Positive-Sine: 
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which, in non-dimensional form, is expressed as 

= -7;*)|sin(m*)| + 7’ / *, 0</* <1 

Exponential Growth / Power Law Decay (m < 0): 
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which is written in non-dimensional form as 


Available Power Minima Under Sawtooth Teet Profile 
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Fig. 4 Gas available power minima curves under 
balanced orbit conditions with a cyclic sawtooth inlet 
temperature profile. 
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For this test case, the system is allowed to reach a 
steady-state mode before decaying into an eclipse. The 
growth constant A is found by matching the initial rate 
of temperature rise from a previous balanced orbit 
mode. For the other two test profiles above, the growth 
constant is calculated to be 

1. Sawtooth: A = 1 


2. Positive-Sine: 


A = 


71 

2 


Figs. 4-6 illustrate the cyclic variation of minimum gas 
available power in response to the cyclic inlet 
temperature profiles outlined above. 


Available Power Minima Under Positive-Sine Test Profile 





Fig. 5 Gas available power minima curves under 
balanced orbit conditions with a cyclic positive-sine 
inlet temperature profile. 
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Dimensionless Gas Inlet Temperature 





which can be regarded as an integral constraint on the 
functions 


T o/T*(t*), and p raln (t*) 


Now consider the maximum SOC given by 


0 = 1 


which corresponds to p : = P max . Therefore, along curves 
of (O.PXI.P™*), the following expression holds: 


B B = 

r max r min 



( 7 
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Fig. 6 Gas available power minima curves under 
balanced orbit conditions with a cyclic exponential 
growth / power law decay inlet temperature profile. 


or equivalently 
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Extrema of Solar Heat Receiver SOC 

First consider the minimum SOC given by 


<D = 0 


which corresponds to p t =l. Therefore, along curves of 
(<D,Pi)=(0,l), the following expression holds: 


P^*)=irQ„ 

Ste 


1 — 




dt 


(u*-s *) 


which is integrated to 


[u * ( t *) -U* (0)]- [s * (t*) -S* (0)] = 





where e is the thermal capacitance ratio 

Me 


e- 


mT cycC p 


In addition, if the system has reached a balanced orbit, 
then the integral of the unsteady term vanishes since 
U*(l) = U*(0) and S*(l) = S*(0), which results in 





where physically meaningful results are obtained when 



>1 


Graphical representations of the maximum SOC are 
shown in Figs. 7-9. 



pmfri 


Fig. 7 Loci for maximum SO C as a function of 
minimum gas available power for selected values of 
receiver available power at the aperture, and a fixed 
combination of thermal capacity parameters. 
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Loci for Maximum 80C (0*1) 



Fig. 8 Loci for maximum SOC as a function of 
minimum gas available power for selected values of 
receiver available power at the aperture, and a fixed but 
higher combination of thermal capacity parameters. 
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Fig. 9 Loci for maximum SOC as a function of a 
combination of thermal capacity parameters for selected 
values of the minimum gas available power, and a fixed 
value for receiver available power at the aperture. 
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Fig. 10 Qualitative illustration of temporal variation of 
incident power crossing the aperture plane. 

Recall that these SOC maxima curves correspond to an 
SD system operating in steady-state mode (SSM). For 
example, this mode can be induced by boosting the 
spacecraft into higher orbital altitudes, which extends 
the sun period and reduces the eclipse period. As might 
be expected, the extended sun period drives the SD 
system in general and the solar heat receiver in 
particular to a state of thermodynamic equilibrium 
wherin the various temperature (and other measurable 
parameters) transients are damped out. Fig. 10 shows 
in a qualitative sense the temporal variation of incident 
power entering the receiver. The discontinuity shown at 
the beginning of each eclipse is not a real effect since 
the actual transition from the sun phase into the eclipse 
phase is a rapid continuous decay, rather than a sharp 
discontinuous drop. However, it is a computationally 
convenient way to model the transition from sun phase 
to eclipse phase. Furthermore, it can be shown that the 
profile shown in Fig. 10 can be generated by a function 
given by 

NCY \ f J ~ l f te 'N 1 

MH Xv, +t ~, Jjj 

where to = 0, NCY is the total number of orbit cycles, 
and H is the Heaviside function. In addition, notice that 
the subscripted parameters allow for variations from 
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cycle to cycle. These variations may be due to the need 
for increased power level, decreased power level, 
increased sun period, or extended eclipse period. In 
addition, the incident power across the aperture may be, 
in general, time-dependent due to time- varying 
shadowing effects on the concentrator or other short 
transients such as concentrator mis-pointing due to 
plume loads from reaction control jets and/or gravity- 
gradient effects. 2 

As pointed out in Hall, III et al. 5 and Mason, 8 the 
anticipated amount of incident power crossing the 
aperture plane is approximately 12.5 kW, and for the 
orbital altitude corresponding to 250 nmi, the total orbit 
period is 93 minutes with about 66 minutes of sun 
exposure and 27 minutes of eclipse. In NASA’s 
Ground Test Demonstration (GTD) system, the solar 
heat receiver uses a eutectic mixture of LiF-CaF 2 as the 
PCM (total mass of 53 lb m or 24.04 kg, heat of fusion of 
340 Btu/lb m or 790 kJ/kg, and melting point of 1873 R 
or 1040 K) and a low-Prandlt-number (for a gas) 
mixture of He/Xe for the working fluid (molecular 
weight of 83.8, c p = 0.059 Btu/lbn/R), the properties of 
which are approximated using ideal gas assumptions. 
The TAC of the Brayton engine is capable of reaching 
speeds of up to 58,000 RPM, with a corresponding 
He/Xe mass flow rate of up to 0.36 lbm/s or 163.3 g/s. 
These numbers correspond to a thermal capacitance 
ratio of £ s 0.2106, Stefan number of Ste = 2.6, and 
dimensionless incident power across the aperture of 
approximately 13.54, assuming a dead state temperature 
of 360 R or 200 K. Also, note that the ratio of Stefan 
number to thermal capacitance ratio is St e/e = 12.3. 
The only other unknown parameter is the virtual source 
or effective cavity temperature, which is a nonlinear 
function of cavity geometric parameters and canister 
surface temperatures. 

Ultimately, these parameters are used to determine 
the maximum SOC corresponding to |} max once the 
turbine requirements are known through the necessary 
minimum gas available power. 


Conclusions 

The theoretical framework for the determination of 
the thermal state -of -charge (SOC) in solar heat 
receivers employing encapsulated phase change storage 
has been developed. The concepts of available power, 
virtual source temperature, and minimum gas available 
power have been used in the underlying theoretical 
analyses. In addition, qualitative and quantitative 
descriptions of minimum and maximum SOC have been 
presented parametrically. Similar parametric curves can 
be generated for non-extremum SOC . 
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